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ABSTRACT 

We present Herschel Space Observatory @ PACS spectra of GQ Lup, a protoplanetary disk in the 
Lupus star-forming region. Through SED fitting from 0.3/im to 1.3 mm, we construct a self-consistent 
model of this system's temperature and density structures, finding that although it is 3 Myr old, its 
dust has not settled to the midplane substantially. The disk has a radial gradient in both the silicate 
dust composition and grain size, with large amorphous grains in the upper layers of the inner disk and 
an enhancement of submicron, crystalline grains in the outer disk. We detect an excess of emission in 
the Herschel PACS B2A band near 63/xm and model it with a combination of ^15 to 70/xm crystalline 
water ice grains with a size distribution consistent with ice recondensation-enhanced grain growth and 
a mass fraction half of that of our solar system. The combination of crystalline water ice and silicates 
in the outer disk is suggestive of disk- wide heating events or planetesimal collisions. If confirmed, this 
would be the first detection of water ice by Herschel. 

Subject headings: Infrared: stars — protoplanetary disks — radiative transfer — stars: individual 
(GQ Lup) 



1. INTRODUCTION 

While they are a small fraction by mass of material in 
the interstellar medium (ISM) and the primordial disks 
around pre-main sequence stars, dust grains play pivotal 
roles in the physical processes that shape these disks. 
The optical pro perties of the dust di ctate the structure 
of the disk (e.g. DAle ssio et al.ll2006l ). and spectral fea- 
tures can be used to probe dynamical processes such 
as settling, turbulence, and grain growth in the disk 
atmosphere ( Bouwman et all 120081 : iWatson et all 120091 : 
iTurner et al.ll2010D . 

In addition to the current disk structure, the grains' 
composition, structure, and size distribution record the 
local thermal history and kinetic conditions of the regions 
in which they formed. For example, silic ate dust grains 
are a t least 98% amorphous in the ISM (jKemper et al.l 
I2004D but are fou nd to be partly crys t alline in pro- 
topla netary disks (|Bouwman et al.l [20031 : ISargent et al.l 
2009). Crystallization is highly sensitive to the specific 
local thermal conditions; different silicate stoichiome- 
tries form at different temperatures and pressures by 
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thermal annealing or su blimation and re-condensation 
(|Henning fc Meeusl 120091 and references therein). Both 
processes occur at temperatures greater than ~ 700K. 
Nominally, this temperature is characteristic of the in- 
ner 1AU of a typical disk surrounding a T Tauri star 
(TTS). However, crystalline silic ates are observed out 
to 40AU in solar system comets dCrovis er et alJll997t ) 
and some TTS systems (jEspaillat et al.l 120071 ) Either 
transportation from the inner 1AU to the ou ter disk via 
winds, turbulent diffusion, or radial flows (jGaill 120011 : 
ICieslal l2009h or in situ formation by heating events 
above the 700K thresh old (shocks or planet destruction) 
(jHarker fe Dcsch 2002) are required to explain these re- 
sults. Characterizing the mineralogical zoning of TTS 
disks at different radii is necessary to discriminate be- 
tween these different scenarios. 

Another major component of dust in disks is wa- 
ter ice, which may play an important role in forming 
planetesimals, eith e r by helping grains stic k together 
i|Ormel et ail 120111: Kuroiwa fc Siron c1 l20Ll or by in- 
creasing the dust to gas mass ratio at the snowline, which 
may induce the accu mulation of grains at this location 
(|Kretke fc Linl 120071 ). Although ice has been detected 
in absorption at NIR wavelengths in highly inclined 
disk systems (|Pontoppidan et all 120051 ). and hot water 
vapor has been det ected in MIR Spitzer spectra (e.g. 
iCarr fc Naiital I2008D , these observations probe mainly 
the upper layers of the disk, leaving the distribution and 
phase of water in the outer disk and midplane largely un- 
known. Recent detections of cold water vapor in TW Hya 
require much less water than indicated by MIR, suggest- 
ing that icy grains may settle to the midplane leaving dr y 
grains in the disk atmosphere (Ho gerheiide et all 120 111) . 
Identifying where and how much ice is present in TTS 
systems would help resolve this question. 

With t his motivation, we investigat e GQ Lup, a 3 
Myr old (Scp eruelo Duarte et al.1 l2008f ) protoplanetary 
K7 TTS system in the Lupus molecular cloud complex 
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with a sub-stellar mas s companion at Of.' 7 separation 
(jNeuhauser et al.l I2005T) a nd a circumprim ary disk that 
appears to be truncated (|Dai et al.l I2010D . Given the 
range o f inclination angles found for the system, ^30° 
to 50° (|Broeg et al.l [20071: ISeperuelo Duarte et all 120081 
hereafter SD08), the separation am ounts to 120 t o 160AU 
at the distance of Lupus, 150 pc (jFrancol [20021 ). Using 
Herschel PACS, we obtained 55 to 145^m spectra of GQ 
Lup. To characterize simultaneously the distribution of 
silicates and water ice in relation to the disk structure, we 
combined these data with archival Spitzer spectroscopy 
and ancillary photometry and used irradiated accretion 
disk models to fit the spectral energy distribution (SED) 
of GQ Lup. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed GQ Lup using Herschel (Pilbratt et al. 
120101 ) on 8 January, 2012 (OBSID 1342238375) with 
PACS (jPoglitsch et aLlfeOlOf ) range spectroscopy modes 
B2A (51-73/zm) and R1S (102-145/um) at Nyquist- 
sampling (R~1500) and a total time of 7774 seconds. 
The data were reduced using t he standar d data reduction 
pipeline in HIPE version 9.0 (|Ottl 120101 ). We extracted 
the spectra from each spaxel, confirmed that the source 
was point-like and well centered on the central spaxel 
within the pointing uncertainty of ^2 , and applied the 
PSF correction to the central spaxel spectrum. The un- 
certainty in PACS absolute flux calibration can be up 
to 30%; however, GQ Lup was observed by both IRAS, 
at 60 and 100/xm, and AKARI, at 65 and 90/jm. We 
use this photometry to confirm the absolute photomet- 
ric accuracy of the PACS spectrum. The point-to-point 
variation of the spectrum after rebinning by a factor of 
10 is ~15%; we assume this as our relative spectral un- 
certainty. 

The Spitzer IRS (lHouck et al.1 l200l low (SL, 
5-14/xm, A/AA=60-120, AORID 5644032) and high 
(SH, 10-19^™, LH, 19-35^m, A/AA=600, AORID 
27064576) spectral resolution data were observed on 30 
August 2004 and 2 September 2008 as part of pro- 
grams 172 and 50641, respectivel y. We reduced them 
wi th SMART (Iffigdon et al.l l200l in the same way as 
m IMcClure et al.l (I2010T) . with the exception that the 
SH/LH data were sky subtracted from off-source frames 
included in that AOR. We estimate the spectrophoto- 
metric uncertainty to be ~5%. 

3. ANALYSIS 

The SED of GQ Lup, is shown in Figure [T] It has 
a strong excess at all infrared wavelengths, indicating 
the presence of a dust sublimation wall and disk. How- 
ever, the disk emission drops off rapidly with increasing 
wavele ngth, consistent with the conclusion by lDai et al.1 
(2010) that it is outwardly truncated. The Herschel B2A 
spectrum shows a peaked triangular shape around 63/im 
suggestive of the water ice feature located there. We 
see no evidence for a forsterite feature at 69^m. In the 
IRS spectrum, we identify the major crystalline features 
by fitting a non-parametric locally weighted scatterplot 
smoothing (LOWESS) baseline to the data, taking this 
as the 'dust continuum' beneath the molecular lines, and 
subtracting a linear fit to regions between known crys- 
talline silicate features to the IRS spectrum (Fig. [2ji). 
There are strong forsterite features at 23 and 33/im, 



blended forsterite-enstatite features at 18 and 28/xm, and 
weak enstatite features around ll^im. 

To determine the composition and structure of the 
disk, we construct temperatur e and density structures 
using the lD'Alessio et all ((2006) irradiated accretion disk 
models, which assume the disk is heated by stellar ir- 
radiation and viscous dissipation. Steady accretion and 
viscosity ar e parametrized through co nstant Mand a, re- 
spectively (jShakura fc Sunva cv 1973). The disk consists 
of gas and dust, the latter of which is comprised of two 
grain populations mixed vertically. Settling is parame- 
terized through e = £/ '^standard, where the denominator 
is the sum of the mass fraction of the different compo- 
nents relative to gas and the numerator is the mass frac- 
tion in the small dust population. 

The silicate and graphite grains have size distributions 
n(a) = noa~ 3 ' 5 , where a is the grain radius with limits of 
0.005/im and a max . To test whether the ice grains have 
grown larger than the silicate grains, we consider three 
size distributions: (Case i) the same power law depen- 
dence and a max as the silicate and graphite grains; (Case 
ii) the same power law dependence but a larger a mQX ; 
(Case Hi) a power law d ependence of n(a) = riQa 2 (e.g. 
iKuroiwa fc Sironoll2011| ) . In this model, ice exists every- 
where below the thermal sublimation temperature; while 
other processes, e.g. photodesorption, are likely to mod- 
ify this snowline (jOberg et al.1 l2009h . a self-consistent 
treatment is beyond the scope of this letter and reserved 
for future work. 

We compute opacities for the gr aphite and water ice 
grai ns using optica l constants from lDraine fc Leel (1984) 
and iWarrenl (|l984) and Mie t heory, assuming th at the 
grains are segreg ated spheres (pollack et al.l ll994). Sili- 
cates are divided into amorphous and crystalline versions 
of two stoichiometries: olivines (Js,lg2~ixF&ixSiO^) and 
pyroxenes (Mgi_ x Fe x SiOz), where x = Fe/{Fe + Mg) 
indicates the iron content. Opacities for the amorphous 
olivi ne and pyroxene are com puted with optical constants 
from lDorschner et al.l (|1995f ) that have x = 0.5 for olivine 
and range from 0.05 to 0.6 for pyroxene. The opacities for 
crystalline olivine (forsterite) and pyroxe ne (enstatite) 
are ta ken directly from those calculated bv lSargent et al.l 
(|2009l) . 

We implement a vertical dust sublimat ion wall with an 
atmos phere following the prescription of lD'Alessio et all 
(2004). Dust is present once the disk temperature drops 
below the dust sublimation temperature, and this radius, 
Rwaii , defines the radial location of the beginning of the 
wall atmosphere, where the dust is optically thin. We al- 
low the dust properties of the wall to vary independently 
of those in the disk to simulate the effects of a radial gra- 
dient in the inner disk mineralogy; in particular, we vary 
the iron content of the silicates in the wall, while in the 
disk we assume the iron content of the bes t-fitting amor- 
phous olivine and pyroxene determined by iSargent et al.l 
( 2009) from comparison with disks in Taurus. 

As input to the code, we assume the stellar and ac- 
cretion properties within the range given by SD08 (val- 
ues listed in Table [IJ . They measured high veiling in 
the system over the course of a month and estimated 
mass accretion rates between 10 -8 and 10 -7 M Q /yr from 
comparing their observed emission lines with magneto- 
spheric accretion models; we adopt an intermediate value 
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of 5xl0 -8 M Q /yr. The outer radius is not totally un- 
constrained, as the mass ratio range of the companion is 
0.014 to 0.03, which is too large for the circumprimary 
disk to radially extend over the 1:2 binary orbital reso- 
nance (|Artvmowicz fc Lubow|[r994l) . Therefore the disk 
can be at most 70AU in radius. 

We create a grid of ^350 models varying e, a, a maXjS , 
the ice mass fraction (f m ,i), silicate composition, and 
Rd and minimize the reduced % 2 statistic, Xri between 
the models and observations to determine the best fit. 
While there are degeneracies between e, a, f m ,i, and Rd, 
the large wavelength coverage of our SED allows us to 
break some of them. For example, only e > 0.1 can 
provide enough continuum emission to fit the IRS spec- 
trum, while e < 0.5 fits the slope between PACS R1S 
and the 1.3 mm photometry. Likewise, smaller disks (e.g. 
Rd ~20AU) underproduce the 20 to 35/im absolute flux 
while larger disks (e.g. Rd ~70AU) overproduce the 100 
to 140/im flux. Once these parameters have been con- 
strained, this implies an a between 0.1 and 0.01. How- 
ever, since the fit to a comes mainly from the submillime- 
ter flux, which depends on the surface density, £ cxM/a, 
the estimated range of Myields 0.02< a <0.2 for the 
same fit to the SED. 

Models without ice fail to reproduce either the shape 
or absolute flux of the PACS B2A spectrum; in Fig [5^ 
we show this model and how the SED changes when ice 
with the distributions given by the previously defined size 
distributions Case and Hi is added. Case i lowers the 
flux too much over the Spitzer range and fails to produce 
enough flux over B2A. Case ii produces only a slight 
increase in the flux between 55 and 70/im. With Case 
Hi, we are able to produce enough flux over both the IRS 
and PACS range. For Case Hi, the B2A and R1S spectra 
are better reproduced by different a maXl i, 15/j,m and 50 
to 70/im, respectively (Fig. [3}d). The mass fraction of 
ice that produces the best fit to the data for this disk 
structure is 2xl0 -3 with respect to gas, or ^25% of the 
total dust content, w hich is half of the solar composition 
(|Pollacket al.l [l99l . Although our model includes ice 
at this abundance in the midplane as well, because the 
disk has so little dust settling, we would require spatially 
resolved submillimctcr photometry to test whether the 
midplane ice abundance differs from that of the upper 
layers. 

The best fitting wall has an a max of 3/mi and is com- 
prised of 100% glassy pyroxene with Fe/(Fe+Mg)=0.4. 
With the assumption of a radially uniform disk com- 
position, the best-fitting models require silicate grains 
with a maa; =0.25/mi and a mixture of 80% amorphous 
silicates, of which 5% are olivine and 75% pyroxene 
(Fe/(Fe+Mg)=0.2), and 20% crystalline silicates, evenly 
split between forstcritc and enstatite. Although our 
model does not yet include radial composition gradients, 
we can still use it to probe the location of the crystalline 
silicates. We compute a grid of theoretical disks trun- 
cated outwardly at 5AU intervals from 50AU to 5AU 
for compositions of pyroxene-forsterite and pyroxene- 
enstatite, allowing the individual crystalline silicate frac- 
tions to vary between 5% and 30%. By isolating the 
model forsterite features in the same manner as for the 
IRS spectrum, we see that the models with sufficient 
forsterite abundance to fit the 33/im feature produce too 



much flux in the 23/im features (Fig. [5]) and vice versa, 
indicating that forsterite is not uniformly distributed 
throughout the disk. 

To compare the data with the emission from an annu- 
lar region of the disk, we subtract the emission of the 
other radii from the 50AU disk and compare the ratio 
of the 23 to 33/im feature to determine which emitting 
annulus is the best match to our data. Forsterite dis- 
tributed between 20 and 50AU fits the relative fluxes 
best and implies a forsterite silicate fraction of ~18% in 
this region. A match to the observed enstatite emission 
of an annulus of 10 to 20AU with a fraction of 15% is ob- 
tained in the same way. We note that the non-detection 
of the 69/tm forsterite feature in the PACS spectrum 
places an upper limit of 90% on the forsterite fraction 
in the outer disk upper layers. To cross-check our re- 
sults, we p erform a two-tempe rature spectral decompo- 
sition (e.g. iSargent et al.ll2009l ) on the IRS spectrum. Its 
results confirm that the inner disk is dominated by large, 
amorphous pyroxene grains with no crystalline grains de- 
tected at > la. The outer disk, however, contains small 
pyroxene grains with forsterite and enstatite present at 
3cr and 2tr, respectively. The resulting system properties, 
assuming radially constant dust abundances in the disk, 
are listed in Table [T] The final best fit has ^=1.2, and 
is displayed in Fig. [TJ 

4. DISCUSSION AND CONCLUSIONS 

4.1. An iron-rich and turbulent inner disk? 

The radius of the wall is proportional to the ratio of 
a grain's opacity weighted by the stellar and local tem- 
perature and the assumed dust sublimation temperature. 
Of the grains we tested, grains with x—QA (40% Fe con- 
tent) had the largest opacity ratio, so they absorb heat 
more efficiently than they emit it and achieve the dust 
sublimation temperature at larger radii than grains with 
less iron. For a wall whose height is given in a fixed num- 
ber of gas pressure scale heights, the solid angle to the 
observer is then larger, and therefore the peak observed 
flux is greater. The high iron content of our inner disk 
is consistent with the Fe/Mg ~1 ratio foun d, on average, 
for p re-solar silicate grains in meteorites (jNguven et al.1 
2010). However, this work does not explore the possi- 
bility that other near-infrared continuum sources that 
might be expected in the inner disk, e.g. metallic iron 
which represents 60% of Mercury's total mass, could be 
present in the wall. The large maximum grain size of 
the wall, taken in conjunction with the high Mand lit- 
tle settling, is indicative of the turbulent enrichment of 
the inner di sk atmosphere with l a rge gr ains seen in other 
systems by ISicilia-Aguilar et al.l (|2007l ). It also suggests 
that there is a radial gradient in the maximum grain size 
of the upper layers, as the outer disk is best fit with 
0.25/im grains. 

4.2. Ice-enhanced grain growth 

Although the silicate and graphite grains are small in 
the outer disk, we can reproduce better the broader shape 
and flux of the SED from 50 to 70/im with a population of 
larger ice grains in the upper disk layers (Fig. (SJd). The 
63/im peak is best fit by an a max of 15/xm and a larger 
e, while a range of a max between 20 and 50/im can pro- 
vide the bulk of the excess flux between 30 and 100/im. 
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The fact that we need ice grains two orders of magni- 
tude larger than the silicate grains and with a different 
grain size distribution power law is consistent with sim- 
ulations for ice enhancemen t of grain growth that t ake 
recondensation into account (jKuroiwa fc Sirondl20ll . if 
we assume that their ice-coated dust grain opacities can 
be approximated by our segregated ice grain opacities 
when a max> i is large. The better fit to the 120/im region 
provided by the 50/im ice grains may indicate either a 
radial or vertical gradient in the ice grain size, which we 
do not yet implement in this version of our model. 

We can think of two potential explanations for why 
we might see an ice emission feature in this system and 
not (to the best of our knowledge) in other systems. 
First, the disk of GQ Lup is highly unsettled (e=0.1- 
0.5). This implies mixing of grain populations between 
the midplane and upper layers. If icy grains are in- 
clined to settle to the midplane or are formed there, we 
may see them because they are being turbulently lifted 
back to the disk surface. Additionally, the lack of set- 
tling means that the location of the Ay — 3 isocon- 
tour (which defines where i ce ceases to be photodesorbed 
(jGorti fc Hollen bach 20Q3)) is at a large enough vertical 
height in the disk that some ice should remain on the 
grains above the r = 1 surface for 63/im emission. Al- 
ternatively, if there are planetesimals in the disk, they 
would not appear to be large enough to clear a gap on 
the order of tens ofAU, as we do not see evidence for one 
in the SED (althou gh we note that this does not con- 
strain smaller gaps; lEspaillat et afll2010t ). The presence 
of a massive, gaseous disk of i?d=30-100AU at the onset 
of the formation of Jovian planets has been simulated 
to scatter smaller planetesimals i n the disk, populat ing 
high inclination orbits at -40 AU (|Kretke et alll2012t ). If 
some of these icy planetesimals were to collide, it could 
generate a reserv oir of icy grains in the upper layers (e.g. 
ILisse et al.ll2012h . 

4.3. Evidence for localized heating in the outer disk? 

The presence of the 63/im feature is also of particular 
interest, as it appears only in crystalline ice, and the 
crystallization temperat ure of initially a morphous water 
ice is typically > 130K (jBaragiolal 120031 ). Although the 
densities in our disk structure imply an ice sublimation 



temperature of at most 120 K (lPollack et al.l 119941 ). it 
has been shown that ice deposited on grains at 160 K and 
subse quently cooled to 1 4 K maintains the 63/im opacity 
peak ([Moore et al.ll200lT ) . This implies that GQ Lup has 
experienced significant heating and cooling events in the 
outer disk. 

Localized heating could also explain the annular crys- 
talline distribution suggested by the last analysis in Sj3j 
This distribution is not indicative of transport models, 
which we would expect to yield monotonically decreasing 
abundances, like those found bv lvan Boekel et al.l (2004) 
for some Herbig AeBe stars. However, there are in situ 
mechanisms that cou ld form crystals at tens of AU. For a 
typical T Tauri star, lHarker fc Descbl (|2002[ ) found that 
shocks from gravitational instabilities can anneal dust 
between 5 and 10AU, assuming an annealing tempera- 
ture around 1200 K. Crystalline dust could also be pro- 
duced thro ugh collisions of p arent bodies, as seen in older 
stars (e.g. ILisse et ail 12012). Given the 3 Myr age of 
GQ Lup, it is not unfeasible that KBO-size planetesi- 
mals could have formed in the outer disk. 

In conclusion, the mineralogical and disk structure 
properties of GQ Lup make it an interesting object with 
which to test theories of planetesimal formation. The 
grains have grown to micron size in the inner disk, while 
the outer disk is highly unsettled and shows signs of water 
ice in its PACS spectrum at sizes consistent with models 
of ice-enhanced grain growth. The presence of crystalline 
ice and silicates at temperatures below their crystalliza- 
tion temperature and in annular distributions is indica- 
tive of a periodic in situ heating mechanism. Follow- 
up studies should be conducted to confirm the disk size, 
resolve any radial structure at larger radii than those 
probed by our spectra, and probe the midplane grain 
size distributions. 
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Fig. 1.— SEP (orange lines and symbols) for GQ Lup. Photometry are from lCovino et al.l JT99^) ■ 2MASS, WISE, AKARI, IRAS, and 
IDai et al.l (1201(1) . Spect ra are from the Spitzer Heritage Archive and this work. The best fitting non-ice model is shown, along with two ice 
models. One fits everything but B2A (50^tm grains, solid grey) and the other fits everything except 120-140/^m (15/mi grains, solid black). 
The remaining model parameters are given in Table [T] The model does not fit the optical data because we do not include emission from 
the accretion shock itself. 
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Fig. 2. — a) Method for isolating crystalline silicate features in the IRS spectrum, b) Observed crystalline features, c) Opacities for three 
most common crystalline silicates, d) Model forsterite features for fixed forsterite mass fraction and varied outer disk radius, as well as 
an annulus. e) Model enstatitc features with best-fitting annulus. f) Best-fitting combination of annuli, 18% forsterite and 15% enstatite. 
The poor match to the 18/xm region is due to fitting a linear baseline to a region that has intrinsic curvature. 
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Fig. 3. — Effects of varying ice grain properties, i.e. mass fraction and size distributions (panel a) and maximum grain size (panel b), for 
a fixed disk structure and silicate/graphite dust properties. 
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TABLE 1 
Stellar and Model Properties 



Parameter 


Value 


Stellar 


(from SD08) 



T ef f 4060 K 

L, l.OL© 

Ay 0.5 mag 

M* O.8M 

M 5 x lO~ 8 Af j/r~ 

i 50° 

d 155±8pc 

Wall 



z wall 

composition 
silicates 
graphite 

Disk 



1500 
0.19AU 
0.03AU 

3.0 

100% pyr, 40% Fe 
0.004 
0.0025 



a 0.01<a<0.1 
e 0.K e <0.5 

R d 50AU 
a max 0.25/im (sil., graph.) 

15 - 70/Ltm (H 2 ice) 

M d O.O12M 
silicate composition: 

(radially constant) 10% forstcrite, 10% enstatite, 

80% pyroxene 

(annular, suggestive) 18% forsterite, 15% enstatite in outer 20-50AU 

abundances 0.004 (sil) 

0.0025 (graf) 
2X10" 3 (H 2 ice) 



